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Abstract 

We study the nucleon axial charge in the chiral perturbation theory for a mixed lattice action of 
Ginsparg- Wilson valence quarks and staggered sea quarks. In particular, we investigate the lattice 
spacing a 2 -dependence of the neutron to proton axial transition matrix element. By using the 
known low-energy constants (LEC's) and an estimated value of a new LEC which appears in the 
calculation, we see a large lattice spacing effect on qa- 
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I. INTRODUCTION 



Lattice QCD, the only known first principle non-perturbative approach to QCD, has made 
a dramatic progress in the last decade because of the increasing power in the computing 
resources and the improvement of the numerical algorithms for the simulations. Today 
fully dynamical lattice QCD simulations together with effective theories enable one to make 
first principle predictions of hadronic physics. Recently, a so-called mixed-action lattice 
simulation using different kind of fermions in the valence and sea sectors has gained attention 
because these simulations can reach the pion masses as low as 300MeV [l|, H, 0, 0, @, @, 



0, @, @, US, ll]- In particular, despite the controversy of the validity of the fourth root 
trick used to overcome the fermion taste doubling problem of the staggered fermions, the 
lattice calculations using Ginsparg- Wilson valence quarks in staggered sea quarks has become 
popular because many quantitative results are obtained [l|, H, UH, @, 0, 0, U|- Further, the 
publicly available MILC configurations 12| have encouraged the use of lattice calculations 
with dynamical staggered fermions. 



The lattice action can be described by Symanzik action [13j, [14j which is based on the 
symmetries of the underlying lattice theory and can be organized in powers of the lattice 
spacing a: 

£ Sym = C + aC^ + a 2 C {6) + ... , (1) 

where represents the contribution from dimension-n operators. The symmetries of the 
lattice action are respected by the Symanzik Lagrangian £sym order-by-order in a. In order 
to address lattice spacing corrections, based on the Symanzik action, several mixed-action 
chiral perturbation theories are developed to investigate the systematic errors which arise 



from lattice simulations due to the non-vanishing lattice spacing [15|, [16j, [17|, [18|, [19|, [20 
These mixed-action chiral perturbation theories provide a way to test the results from the 
lattice calculations using the mixed lattice action and vice-versa. 

The nucleon axial charge, gA, is an important quantity in QCD which quantifies the 
spontaneous chiral symmetry breaking. Since its value is known to a very high precision in 
experiments 
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it can be used as a very good test for the first principle lattice calculations. 
This quantity of hadronic physics has been studied extensively in both lattice simulations and 
chiral perturbation theory [!, 22, 23, 24, 25, 26, 27|, H, 2§]. Although to perform the lattice 
simulations to calculate g\ is straightforward, some controversy remains. For example, the 
gA obtained from lattice simulations has raised the controversy of the possibility of large 
volume effects in the chiral limit [30|, 31]. Recently, using the mixed-action of Domain Wall 
valence quarks with staggered sea quarks, [5J obtained a nice result of g& which seems to 
match nicely with the prediction from chiral perturbation theory and it is expected that 
the future data might be able to make contact with experiment via extrapolation. On the 
other hand, in [5|, the discretization errors in computing g& using the mixed-action might 
be overlooked. With the advent of mixed-action ;^PT [19j and extension to the baryon 
sector |2C|, it is surprising that lattice spacing effects on the nucleon axial charge have not 
been addressed. In this paper we address this issue mentioned above, namely, we use the 
mixed-action chiral perturbation theory to study the lattice-spacing dependence of In 
particular, we found that it depends on the unphysical low energy constant g\ and a new 
unknown LEC for which the actual value must be obtained from lattice calculations. 

This paper is organized as follow. In Section 2, we briefly review the chiral perturbation 
theory of Ginsparg- Wilson valence quarks in staggered sea. In particular, we will put em- 
phasis on the lattice-spacing a dependence of the mixed- action Lagrangian. Next, in Section 
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3, we will construct the axial-vector current from the mixed-action Lagrangian and compute 
the pn matrix element of the axial- vector current. Further we will compare the results ob- 
tained in the presence and absence of lattice-spacing a. Finally we conclude in Section 4. 
For completeness, we include some standard functions which arise in this calculation in the 
Appendix. 



II. CHIRAL LAGRANGIAN 



To address the finite lattice spacing issue, with the pioneering work in 32, 33], 15, 16|, 34 
have extended the in the meson sector. This finite lattice spacing artifacts has also 
be investigated in staggered xPT for the mesons 35j, |3g, |37j and the mixed-action PQ;\PT 
13, H, 19, 20]. To address the finite lattice spacing effects, one utilizes a dual expansion in 
the quark masses and lattice spacing with the usual energy scale [l6|, 17]: 



m, 



and the following power counting scheme [2 



e 2 ~ 



mJA 



QCD 



n 2 \ 2 
a yl QC*D 



(2) 



(3) 



which is relevant for current improved staggered quarks simulations. 

In this section, we briefly introduce the chiral perturbation theory of the partially- 
quenched mixed-action theory. In particular, we focus on the mixed-action of Ginsparg- 
Wilson valence quarks in staggered sea quarks. We only write the chiral Lagrangian of the 
associated mixed-action theory and do not go into the detail of how this chiral Lagrangian 
can be constructed from the Symanzik Lagrangian since the detail of the procedure can be 
found in the references cited above. For our purpose in this paper, we only need to keep in 
mind that = and the taste-symmetry breaking and 50(4) breaking operators will not 
enter at the order we calculate the axial- vector current matrix element 1 201 . 



A. Mesons 

In the following, the strange quark mass is assumed to be fixed to its physical value, 
therefore one can use two-flavor theory and does not need to worry about the extrapolation 
in the strange quark mass. The lattice action we consider here is built from 2 flavors of 
Ginsparg- Wilson valence quarks and 2 flavors of staggered sea quarks. In the continuum 
limit, the Lagrangian is just the partially quenched Lagrangian which is given by: 

£ = QIpQ + Qm q Q , (4) 

where the quark fields appears as a vector Q with entries given by: 

Q = (u,d,j 1 ,j 2 ,j3,ji,h,h,h,h,u,d,) T . (5) 

and transforms in the fundamental representation of the graded group S77(10|2). Notice the 
fermion doubling has produced four tastes for each flavor (j, /) of staggered sea quark. The 
partially quenched generalization of the mass matrix m q in the isospin limit is given by: 

m q = diag(m u , m u , m,£/, rrijCi, m u , m u ) , (6) 
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with as the 4x4 taste identity matrix. The low-energy effective theory of the theory we 
consider above is written in terms of the pseudo-Goldstone mesons emerging from sponta- 
neous symmetry breaking which are realized non-linearly in an C/(10|2) matrix exponential 

E .exp(^)^. (7) 

To the order (next-to-leading order) we work in investigating the a 2 -dependence of the matrix 
element of the axial- vector current between nucleons, the relevant partially-quenched chiral 
perturbation theory Lagrangian for the mesons up to order 0(e 2 ) is given by: 

f 2 1 
£ = ^ str (3^0^) - X str (m 9 S t + mjs) + (str $) 2 + a 2 V . (8) 



where 



/ = 132 MeV, the str() denotes a graded flavor trace and E is defined in (j7]). The M, M, and 
Y are matrices of pseudo-Goldstone bosons and pseudo-Goldstone fermions, for example, see 
The potential V contains the effects of dimension-6 operators in the Symanzik action 



19(. Expanding the Lagrangian in Eq. (jSJ) to the leading order, one can determine the 
meson masses needed for the calculations of baryon observables. In particular, the relevant 
mesons needed for the axial- vector current matrix element calcualtions are the valence pion, 
mesons made of S{V with one staggered sea quark Si of flavor S and quark taste % and a 
Ginsparg- Wilson valence quark V and finally mesons with two staggered quarks in a flavor- 
neutral, taste-singlet combination. The associated masses to the lowest order for the latter 
two mesons can be written in terms of the valence pion mass m vv which can be determined 
from the valence spectroscopy and the pseudoscalar taste pion mass m 2 ss which the mass 
can be learned from the MILC spectroscopy and are given by: 



m 



m 2 ss + rriy V 16 a 2 C n 



sv — 2 f 2 
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™ 2 ss,i = ™ 2 ss + ^r(c 3 + c 4 ) , (10) 

where C m i x and C3 + C4 are the parameters in the potential V. Since these masses are 
independent of the quark taste, we do not specify the taste index in ffTUj) . Notice the flavor 
singlet field $ is rendered heavy by the U(1)a anomaly and can been integrated out in 
PQxPT. However, the propagator of the flavor-neutral field deviates from a simple pole 
form [3^] . Since only the valence- valence flavor- neutral propagators are needed for our later 
calculations, for V, V = u, d, the leading-order r}yr]v> propagator in the isospin limit is given 

iSyy, % (q 2 -m 2 SJ ) 

Q - ™>vv + ie 1 [q - myy + le) 

One can further show that these propagators can be conveniently rewritten in a compact 
form which will be useful in our later calculations: 

^vvn V ' = bvv'Pv + ^ivv {Pvi Pv) ) (12) 
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where 



v 



q 2 - my V + ie 



H VV (A,A) = - 



In writing f|T3|) . we have used A = A{m 



vv) 



SSJ 



m 



vv, 



d 

drn vv 



A - A 



(13) 



B. Baryon 



As has been shown in 20j, to 0(e 2 ), the free Lagrangian for the 572-dimensional super- 
multiplet and the 340-dimensional super-mult iplet fields in the mixed-action 

SU(10\2) partially quenched has the same form as in quenched and partially quenched 
theories 0, [H, Ej with the addition of new lattice-spacing dependent terms: 



C 



(Bv ■ VB) + 2afj Q) (BBM+) + 2(3™ (BM+B) + 2afj Q) (BB) str (M+) + a 2 V B 

(T% ■ VTJ + A (7%) - 2^ Q) (T»M + %) - 2^ Q) (T%) str (M + ) + a 2 V r . 

(14) 



The baryon potentials Vg and Vr in fTL4")) arise from the operators in of the Symanzik 
Lagrangian. In the baryon Lagrangian, the mass operator is defined by: 



A4+ = -(eW f +ZrnQ£). 



(15) 



and the parameter A ~ eA x is the mass splitting between the 572 and 340 in the chiral 
limit. The parenthesis notation used in Eq. (fl4")) is that of j4o[. Further, the embeding 



of the octet and decuplet baryons in their super-multiplets is the same as before [28|, [39 
The Lagrangian describing the interactions of the B^ k and T^ k with the pseudo-Goldstone 



mesons in the mixed-action is again the same as in the partially quenched theories [28] with 
a new a 2 -dependence term: 



C = 2a (BS^BA^) + 2(3 (BS^A^B) + 2H (T^S^A^) 
+ ^C[(TA„B) + (BA u T»)]+a^ t . 



(16) 



The axial- vector and vector meson fields A^ and are defined by: A^ = | (£<9 M £t — 

and V = \ (£(9^ + ft^f). The latter appears in Eq. (TH)) for the covariant derivatives of 

Bijk and %j k that both have the form 



{V»B) ijk = d»B ijk + {V»)% jk + 



<l(Vi+Vj)(Vk+Vn) ( ) ™ Bi n 



(17) 



The vector is the covariant spin operator [221 . |23| and a 2 C\ n \ is from the operators in 



£( 6 ) of the Symanzik Lagrangian. As we will see later, the explicit form of this term is not 
required in our calculations. The effective axial-vector current from a 2 C,f^ t can be obtained 
by a simple argument. The parameters that appear in the mixed-action PQxPT Lagrangian 
can be related to those in x?T by matching. Further, since QCD is contained in the fourth- 
root of the sea-sector of the theory, one should restrict oneself to one taste for each flavor 
of staggered sea quark when performing the matching. To be more specific, one restricts 
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to the qsQsQs sector [20j and compares the mixed-action PQ%PT Lagrangian obtained with 



the xPT. With this identification and matching procedure, one finds [28 



(PQ) 

M 



(PQ) 



(PQ) 

M 



+ Q a M 



(PQ) 



+ - 3 ih, 



Further, when restricting to the tree level, from ffTB]) one also finds 28| 



(PQ) 



(PQ) 

M 



a 



4 1 
- 3 9A + -9! 



ft 



3 91 ~ 3 9A 



n 



9aa , C — 9 AN 



(18) 



(19) 



and gx = 0. 



III. THE AXIAL- VECTOR CURRENT 

The matrix element of the axial- vector current, qT a 'y l ,.y ^q, have been studied extensively 



both on the lattice @, 0, Q and %PT @, H 0, O, [3, For our convention 



we 

will follow [4JJ, |42| and use the following charge matrix for the flavor- changing current in 
extending the isovector axial current Qf a ^^^Q to PQQCD since we are interested in the 
neutron to proton axial transition: 

r + = (0,l,0,0,0,....,0). (20) 

With this convention, at leading order, the flavor- changing axial current is given as [281 ] : 



(PQ) 



Ju,5 



2a (BS^B r+ ) + 2(3 (BS^ rt + B) + 2H (T^ r+T v ) 



\ + f V+f ) and is obtained from £) n ' t m 



where — 2 
insertion of the mass matrix will be at 0(e 4 



a Jafj.,5 



,(6) 



(21) 



161) . Notice since the 
we only need to take the tree level contribution 
from a 2 J^ 5 into our calculation. Therefore, although there are a lot of terms in the axial 
current arise from the a 2 corrections to the lattice operator itself, the net contributions to 
the current 5 can be effectively written as: 



jafj.,5 ~ ^Ca 



(BS,B r+ ) + 2C' a (BS, r+B) 



(22) 



The calculation of the matrix element of the neutron to proton axial transition at next- 
to- leading order is 0(e 3 ) in our power counting. These leading non-analytic terms are from 
the one-loop diagrams shown in Fig[TJ To obtain the complete 0(e 3 ) calculation, from our 
power counting scheme: e 2 ~ a 2 AqCD and e 2 ~ m q / 'Aqcd, we see that in addition to taking 
the a 2 -dependence of the loop meson masses into our calculation, we must evaluate the 5 
at tree level. After carrying out the calculation, one finds: 



pn 



2U pSJJn , 



(23) 
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N ■ .' N 



N N N N 



FIG. 1: The one- loop diagrams which contribute to the leading non- analytic terms of the nucleon 
axial transition matrix element. Mesons are represented by a dashed line while the single and 
double lines are the symbols for a nucleon and a decuplet respectively. The solid circle is an 
insertion of the axial current operator and the solid squares are the couplings given in (|16j) . The 
wave function renormalization diagrams are depicted by the two diagrams at the bottom row. 

where from the contributions of local counterterms involving one insertion of mass 

matrix wq 1 and T pn is given by: 



- g 2 A gi(2Ri(m vv , fx) - 2R 1 (m vs ,n) - 6A(m S s,i)j 

- 9Agl^—Rx{mw^) - Y R i( m vs,n) - 3A(m ss ,i 
4 / 

- 9AN9ig yN^mw, A,/i) - Ni(m vs ,A,n) 



+ gl N 9AA— (ji(m vv , A,/i) + -Ji{m vs , A, fi) 
~ 9an9a ^ -^-^Ni(myy, A, fi) + N\{rays, A, /i) 

- 2^J 1 (m vv ,A,fj l ) + Ji(m yS) A,//)| 



(24) 



1 These local terms have the expressions aim vv + a2fn 2 ss with ii and 0,2 must being determined from the 
lattice calculations 
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FIG. 2: T pn — * qa in the chiral limit. Here the physical value of qa is set to be 1.25. The band 
structure indicates the variation of qa due to the bounds on g^A and gAN- 



where C = C a + C' a , the functions Ji(m, A, /j,), Rx(m, //), Ni(m, A, //) are denned in the 
Appendix, A is the A-nucleon mass splitting, m a b are given in (JTU1) and lastly, the function 
A{ m ss,l) is given by: 



Mmssj) = 3^ (m 2 SS)J - m^) [ log ^ + 1 ] (25) 



2 

l> 2 

All of the couplings in eq. ff24|) take their chiral-limit values. It is easy to see that all 
the a 2 -dependence in f f24|) is contained in m^s, m ss,i and a 2 C. Notice with our Ri(m,fi), 
Ji(m, A,/i) and Nx(m, A, fi) we have subtracted off the chiral and continuum limit values 
of the loop diagrams by hand. This corresponds to a renormalization of the tree level 
coefficients, and produces Qa which is the chiral limit value (Fig [2]). 

The lattice spacing dependence in Eq. ( 12^1) is completely determined by three parame- 
ters, namely, C m i x and C3 + C4, and a new unknown low energy constant C . In order to 
investigate the a 2 -dependence of np axial transition matrix element, in Fig [31 we have taken 
the unquenched limit in which myy = mss 2 and plotted (|24|) at two different lattice spacing 
a = 0.12fm and a = with the pion mass varying from lMeV to 500MeV 3 . In the figure, 
gA is set to be 1.25 and the low-energy constants gi, g<\N, 3aa are allowed to vary within 
their reasonably known bounds:: — 1 < g\ < 1, 1.0 < I^atvI < 2.0 and 2.5 < |#aa| < 3.5. 



Further, C mix is estimated to be 4.5fm < C mix < 9.5fm 431 and G% + Ca is determined 



as C3 + C4 = 2.34fm 6 [l|. Finally, since C is assumed to be of natural size which is 



2 The unphysical effects of the mixed-action PQ theory arise from two sources. One is from the mass 
difference between valence and sea quarks and the other is because of the non-vanishing lattice spacing. 
By taking the unquenched limit, namely, m V y = rn? SSl one can eliminate part of the unphysical effects. 
Notice when the lattice spacing is zero, the physics is recovered from the PQ theory in the unquenched 
limit. 

3 We have dropped the local terms c pn 's since the loop contributions formally dominate over these terms 
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100 200 300 400 500 

(MeV) 

FIG. 3: In this figure, we have plotted T pn as a function of m^. The green band stands for the 
band with a = while the blue band are the T pn obtained from qa = 1.25, C3 + C4 = 2.34fm -6 , 
reasonably known bounds on gi, g&N and g^A- — 1 < 9i < 1, 1-0 < |<?Aiv| — 2.0, 2.5 < |#aa| < 3.5, 
a natural estimated value for C m \ x : 4.5fm~ 6 < C m - lx < 9.5fm~ 6 , and finally a natural variation of 
C'\ — 3fm~ 2 <C'< 3fm~ 2 . The lattice spacing a for the blue band is fixed to be 0.12fm. 

Aq CD ~ 2.3fm~ 2 , we have used — 3fm -2 < C < 3fm -2 for the figure. However, keep in mind 
that the actual value of C must be determined from lattice calculations. 

From the figure, we indeed see a large lattice spacing dependence for T pn . One might be 
concerned that the correction at the chiral limit is sufficiently large that the prediction 
might be breaking down. Further, the band at the chiral limit does not cover the value of 1.25 
which is the expected qa at m n ~ 0. We point out that with a = 0.12fm, the mass square 
corrections to the VS and SS, I mesons are around (400MeV) 2 to (450MeV) 2 which implies 
t licit clt m n ~ on the figure is effectively similar to that obtained at ~ 450MeV 
with a = 0. Therefore a large correction is no suprise. By reducing the lattice spacing a to 
0.08fm, which might be the standard lattice spacing in future simulations, the correction is 
around (270MeV) 2 . Notice by shifting the a = result to the left by 450MeV units, now 
the green band indeed overlaps with the blue band at ~ 0. 

In addition to the large a 2 mass shifts, the lack of unitarity in the PQ theory will also lead 
to divergences when approaching the chiral limit and hence contributes a large correction 
to at m n ~ 0. It is clear that the A(mss,i)s' in Eq. fl24"j) are responsible for this 
divergence since they are from the double pole of the flavor neutral propagator which violate 
the unitarity. Indeed if we plot Eq. (jUj) as a function of m n with C mix = 0, we see a large 
correction to T pn at m v ~ (FigHJ). A closer look at Eq. (|24|) and Eq. f[2"5"j) shows that the 
divergence arises from the difference of valence and taste-singlet pion masses. In principle, 
myy can be tuned to m SS j [3] such that the divergence is less severe. However the large 
m ss,i will still cast doubt on the chiral expansion due to the large mass. Also notice if 
Cmix = C3 + C4 = in Eq. ([21]), then a narrow band centered around the result obtained by 
setting a = should be observed. Indeed this is comfirmed in Fig [5j 

Lastly, to investigate quantitatively the lattice spacing effects on g^, let us focus on 
~ 320MeV which is the relevant smallest pion mass used in most recent mixed-action 
simulations. In Fig [3} we see with a = 0.12fm, T pn receives a 62 percent correction at 
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500 



FIG. 4: T pn obtained by letting C mix = in fl24"j) . 



~ 320MeV. This observation shows that the lattice spacing effects on g& of current 
related simulations should not be overlooked and that it is not clear whether ;^PT can be 
reliably used to account precisely for lattice spacing effects. However, with lattice spacing 
a = 0.08fm, we observe a reasonable 28 percent correction to Y pn at ~ 320MeV. We 
also point out that there is an unphysical parameter gi-dependence in ( )24l) even in the 
unquenched limit. This is a " partially quenched" artifact depending on the mixed- action 
and clearly from (T2"4~l) the gi-dependence will disappear when the lattice spacing a is set to 
zero. 



IV. CONCLUSION 

The mixed-actions provide us with a powerful tool in lattice simulations because they 
enable one to use different type of fermions in the valence and sea sectors. In particular, 
current lattice calcualtions using Ginsparg- Wilson valence quarks in staggered sea quarks 
can reach the dynamical pion masses as low as 300MeV. 

In this paper, we have calculated the pn matrix element of the axial-vector current up 
to 0(e 3 ) order using the mixed-action of Ginsparg- Wilson valence quarks in staggered sea 
quarks. Further, we have detailed the lattice spacing artifacts for this matrix element. To 
0(e 3 ), we found that the pn axial-current matrix element depends on the lattice spacing 
via three parameters, namely, C m ; x , C3 + C4 and finally a new low energy constant C . The 
low energy constant C m { x affects the masses of mesons made from a Ginsparg- Wilson quark 
and a staggered quark. Its physical value can be fixed either from mixed meson masses 



or the pion charge radius 43]. Further, the combination of parameters C3 + C4 has been 
already constrained from staggered meson lattice data [lj]. The new low energy constant 
C, on the other hand, can be evaluated from the lattice simulations of determining the 
nucleon axial charge ja, 0, 25[. However, as has been already demonstrated in SI, the 
continuum extrapolation with only one lattice spacing available will lead to a large amount 
of uncertainty in the physical values of the associated low energy constants. Ideally, this low 
energy constant C will be more accurately determined if a variety of lattice spacings in the 
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FIG. 5: The blue band is the T pn obtained by letting C m i x = C3 + C4 = 0. The green band 
represents the T pn with o = 0. 

lattice simulations are available. 

Lastly, since the finite volume effects are small at the scale of dynamical quark masses 
and the box size available in today's simulations as indicated in jEl, [2^], the formulas given 
here should be sufficient for the comparisons between the predictions from and the 
results from numerical simulations. 
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V. APPENDIX 

In this Appendix, we list the functions needed in our calculations: 



1 



m 2 log 




(26) 
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KAm, A, n) = -I^-A 2 logf^) +-^( (m 2 - -A 2 ) log 
v ' ,M 2 16vr 2 & V4AV 4167T 2 IV 3 / & V /i 2 



3 VA + v / A 2 -m 2 + ze/ 



2m 2 / r-^ - / A - y/A 2 - m 2 + ie\ \ i 

+ — - Trm - VA 2 - m 2 log = J J . 27 

3 A V VA + v / A 2 -m 2 + ze/ ' 1 



Mm.A.u) = --—^A 2 \og(-^—) + ~-^—( 

u ' ' w 216tt 2 & V4A 2 / 416vr 2 l 



m 2 -2A 2 Uog 



+ 2AyA 2 ^log f AWA2 - m2 + ") ] . (28) 
VA + v / A 2 -m 2 + i e/ 1 
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